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摘要 : 重 离 子 碰撞 过 程 中 轻 粒 子 的 产生 对 于 提取 核 物 质 状态 方程 的 
子 分 子 动力 学 (UrQMD) 模型 ， 利 用 统计 衰变 模型 
区 Au+Anu 碰撞 过 程 中 轻 粒 子 的 集体 流 和 核 阻止 本 领 的 影 
该 的 部 分 动力 学 性 质 ， 在 考虑 
量 的 影响 随 碰 撞 能 量 的 升 高 而 减弱 。 结 果 表 明 
子 的 产生 对 敏感 于 核 物 质 状 态 方程 的 观测 量 有 着 一 定 的 影响 ， 在 夭 
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信息 ， 例 如 自由 核子 和 碎片 的 动量 分 布 ， 
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学 演化 过 程 。 一 般 来 说 我 们 可 以 借助 输 运 模型 模拟 重 离 
子 碰撞 ， 并 将 模拟 结果 与 实验 数据 对 比 ， 从 而 提取 核 物 
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信息 有 着 重要 的 作用 。 基 于 极端 相对 论 量 
GEMINI++ 处 理 预 激发 碎片 的 衰变 ， 研究 了 预 碎片 的 退 
向 。 发现 由 于 记忆 效应 ， 重 离子 碰撞 过 程 中 
预 碎片 的 衰变 后 可 以 更 好 地 描述 实验 数据 ， 
重 离子 碰撞 过 程 中 预 碎片 的 衰变 以 及 轻 粒 
用 这 些 观测 量 提取 核 物 质 状 态 方 程 的 信息 
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， 输 运 模型 结合 聚合 模型 得 到 的 预 碎 片 通常 处 于 激 
发 态 ， 而 真实 重 离子 碰撞 过 程 中 碎 
程 中 会 经 历时 间 尺 度 较 长 的 退 激发 。 


7 析 轻 粒子 的 一 些 


片 在 到 达 探 测 器 的 过 
因此 ， 在 相对 较 低 
片 的 激发 能 与 束 流 
能 量 相当 ， 我 们 需要 重视 这 种 预 碎片 的 退 激 对 最 终 的 观 


男 一 方面 ， 对 轻 碎 片 产 额 的 低估 是 输 运 模 型 模拟 中 


能 重 离子 碰撞 中 一 个 长 期 存在 的 问题 史 9, 理论 上 已 经 提 
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9 由 质子 和 和 气 核 的 集体 流 的 结果 。 
为 了 更 好 地 给 出 轻 粒 子 的 集体 流 分 布 以 及 探究 衰变 对 中 
能 重 离子 磁 撞 过 程 中 轻 粒 子 观 测量 
计 模 型 GEMINI++ 月 
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中 能 重 离子 碰撞 ， 并 使 用 统计 模型 一 -GEMINI++ 来 处 
里 预 激发 碎片 的 衰变 ， 研 究 预 碎片 的 衰变 对 轻 粒 子 产 额 Cag ~ FP(p, p) x as， (5) 
4h. 集体 流 和 核 阻 止 的 影响 。 文 章 结构 安排 如 下 : 在 
第 2 节 中 对 使 用 的 理论 模型 和 观测 量 进行 曾 述 ， 第 3 75 
E Pyy > 1GeV/c 
为 结果 与 讨论 ， 最 后 在 第 4 节 给 出 总 结 。 


2 理论 模型 及 观测 量 
动力 学 模型 : UrQMD 


UrQMD 模型 是 一 种 处 理 微观 非 平衡 动力 学 过 程 的 
输 运 模型 821， 在 模型 中 每 个 核子 由 一 个 具有 确定 宽度 


2.1 


的 高 斯 型 ! 波 包 表示 。 在 Au+Au 磁 撞 中 ， 波 包 宽度 参数 
Mor = 2fm， 每 个 核子 的 坐标 和 动量 的 时 间 演 化 遵循 


哈密 顿 正则 方程 。 系统 的 总 哈密 顿 量 (A) 包括 动能 TT 和 
核子 -核子 有 效 相互 作用 势 U， 其 中 U 又 包含 Skyrme # 
能 项 以、 库伦 能 Ucou， 以 及 动量 依赖 势能 项 Ung 
= T + U, + U coul + U mdo 
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撞 ， 本 文采 用 改进 的 量 
中 所 使 用 的 Skyrme 能 量 
平均 场 。 局 域 项 和 动量 相 
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其 中 65 = (Pa 一 Pp) / (Pa + Pp) 为 同位 旋 不 对 称 度 ， pp 为 
质子 密度 ，p, 为 中 子 密度 。 在 本 工作 中 ， 核 物质 的 饱和 
密度 为 pn = 0.16 fm 飞 ， 势 能 参数 分 别 为 a = —393 MeV, 
p = 320MeV,y =1.14,gs = 19.5 MeV fm? ,gsiso = 一 11.3 
MeV fm?, A... =20.4 MeV, Bsn = 10.8 MeV,C,,,, = —9.3 
MeV, n=1.3, tng = 1.57 MeV, apa = 500 (GeV/c) 2, 和 
c = 0。 这 些 参 数 设 置 对 应 一 个 软 的 动量 依赖 的 状态 方程 
(SM-EoS), 相应 的 核 物 质 不 可 压缩 系数 为 Ko = 200 MeV 
、 对 称 能 斜率 参数 为 L = 80.95 MeV. 

对 于 核子 -核子 弹性 截面 ， 我 们 对 其 考虑 了 介质 修 
正 效应 。 通 过 一 个 密度 和 动量 依赖 的 参数 化 修正 因子 
F(p, p) 与 自由 核子 -核子 截面 Of 之 积 得 到 介质 内 的 核 
子 -核子 弹性 截面 64。 
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此 外 ,由 于 核子 为 费 米子 , 满足 泡 利 不 相 容 原理 , 因 
此 在 碰撞 过 程 必须 考虑 Pauli 阻塞 效应 [2 人 I。 散射 过 程 
中 ， 未 态 相 空间 如 果 没 有 被 占据 ， 则 该 散射 为 成 功 的 散 
射 ， 成 功 的 概率 是 (1-f), 否则 该 散射 被 阻 止 。 核子 -核子 
碰撞 过 程 中 出 射 粒 子 的 相 空 间 密 度 分 别 为 f; 和 f: 
(wm) Quim) 20 
fii = : 2, ca m a ; (7) 
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作用 处 于 同一 个 碎片 。 故 而 质 
于 其 他 核子 对 的 值 ， 
虽然 这 样 的 改进 不 能 完全 解决 
HEIRS, 
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论 的 统计 误 


片 退 激 过 程 中 发 生 
的 衰变 ， 包 括 轻 粒 子 的 蒸发 、 对 称 和 不 对 称 裂 变 以 及 其 


型 已 被 广泛 的 用 
撞 过 程 中 的 预 碎 片 的 
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的 参数 为 默认 
壳 平滑 能 级 密度 参数 为 : ky = 7.3 MeV, 
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模型 中 包含 4 个 基本 的 输入 量 : 激发 核 的 质量 数 A, 
EMA Z， 激 发 能 E*， 以 及 角 动 量 L。 其 中 激发 能 可 以 
由 预 碎片 的 结合 能 和 实验 上 测量 得 到 的 基态 结合 能 相 减 


得 到 : 


* _ mexcited ground 
E = Exina T Evind ] (10) 


其 中 ES"! 可 以 从 国际 原子 质量 评估 AME2020 数据 
W68 中 得 到 ， 对 于 数据 表 之 外 的 未 知 核 素 可 以 利用 
液 滴 模 型 (9 计算 其 基态 结合 能 。 当 激发 能 小 于 等 于 0 
时 ， 预 碎片 将 被 视 为 处 于 基态 ， 不 对 其 考虑 衰变 。 预 碎 
片 的 角 动 量 则 可 以 通过 经 典 力学 得 到 [al; 


L= xP, (11) 
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r, 和 p, 分 别 是 预 碎片 中 第 ;个 核子 的 坐标 和 动量 矢量 。 
经 验证 ， 模 型 中 是 否 考虑 角 动 量 的 贡献 并 不 会 对 结果 产 
生 明显 的 影响 。 
2.3 观测 量 

在 本 文中 ， 我 们 主要 讨论 预 碎片 的 衰变 对 经 粒子 的 
集体 流 和 核 阻止 本 领 的 影响 。 集 体 流 中 常用 的 直接 流 0, 
和 椭圆 流 v 分别 定义 为 al: 


2 2 
anf zs ) «(5 2). (12) 
Pi + PS Heres 
其 中 p, All p, 分 别 是 粒子 沿 x 方向 HRSS My A 
向 的 动量 ，p, 则 是 沿 束 流 方向 z 的 动量 ， 同 时 可 以 得 到 
横 动量 p, = /成 + 成， 尖 括 号 表示 对 所 有 事件 中 观测 六 
子 求 平均 值 。 核 阻止 本 领 通常 使 用 INDRA 实验 组 提出 
的 R, I"! 和 FOPI 实验 组 提出 的 varti 0? 表征 : 


2 
ET vartl — o) (13) 


QR» 
EH, E, (E) 是 质心 系 内 横向 纵向 ) HER, RA 


Au+Au Ew = 150 MeV/nucleon b = 0-2 fm 


10 上 * INDRA 
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charge Z 
到 1 (EREE) 束 流 能 量 为 150 MeV/nucleon 时 Au+Au 中 


心 碰撞 (b = 0-2 fm) 中 200 fm/c 时 刻下 考虑 衰变 前 后 
的 电荷 多 重 数 分 布 并 与 INDRA 实验 值 26 对 比 


与 INDRA 实验 值 进行 了 对 比 ， 这 里 采用 的 演化 停止 时 
间 为 200 fm/c。 图 中 蓝 色 空心 圆圈 为 不 考虑 衰变 的 结果 ， 
绿色 实心 三 角 为 考虑 衰变 后 的 结果 ， 红 色 五 角 星 为 实验 
测量 值 。 可 以 看 出 ， 在 考虑 激发 预 碎片 的 衰变 后 ， 由 于 
重 质量 碎片 的 衰变 ，Z > 2 的 碎片 的 产 额 被 整体 的 压低 
并 且 更 加 靠近 实验 值 ， 同 时 Z = 1 和 Z = 2 粒子 的 产 额 
得 到 了 提升 并 且 可 以 较 好 地 描述 Z = 2 粒子 产 额 的 实验 
数据 。 

进一步 的 , 图 2 展示 了 入 射 能 量 为 150 MeV/nucleon, 
磁 撞 参数 为 = 0-2 fm 的 Au+Anu fife ee SCR] f 38 
(p, d, t) MARME GHe 和 4He) 产 额 的 快 度 分 布 , 其 中 
蓝 色 虚线 为 不 考虑 衰变 的 结果 ， 绿 色 实 线 为 考虑 衰变 后 
的 结果 , 左 中 右 部 分 分 别 表 示 在 1= 200, 250 和 300 fm/c 
时 是 否 考虑 衰变 后 的 结果 。 当 考虑 预 激 发 碎片 的 次 级 误 
变 后 p, d 的 产 额 稍 有 增加 ， 氢 同位 素 的 产 额 增益 尤其 明 


盖 所 有 考虑 的 粒子 。 粒 子 在 x 和 z 方向 分 别 具 有 快 度 y、 
和 y,， 其 对 应 的 方差 分 别 为 (六) 和 G2). 


E 
a (14) 
B 2 — Pxz 
入 射 弹 核 的 快 度 大 小 为 wwe = ine, GE SORTE x 


和 z 方 向 的 约 化 快 度 为 yo = yy Yom VelVoror 以 上 
均 指 质心 系 。 

3 结果 与 讨论 

31 产 额 分 布 


首先 ， 图 1 给 出 了 束 流 能 量 为 150 MeV/nucleon 时 
的 Au+Au 中 心 页 撞 考 虑 衰变 前 后 的 电荷 多 重 数 分 布 ， 


显 。 然 而 对 1 的 产 额 快 度 分 布 几 乎 没有 影响 ， 这 是 因为 
的 比 结合 能 较 大 ， 本 身 不 容易 发 生 衰变 ， 以 及 大 质量 的 
预 激 发 碎片 在 蒸发 相同 质量 轻 核 时 相 比 于 发 射 ! 更 倾向 
于 发 射 比 结合 能 较 小 的 3He。 由 图 1 可 知 , 由 于 发 生 衰变 ， 
大 质量 碎片 将 衰变 为 轻 质量 碎片 ， 导 致 轻 粒 子 产 额 分 布 
的 提升 。 因 此 ， 考 虑 衰变 后 He 和 “He 的 产 额 分 别提 升 
为 之 前 的 1.4 倍 和 4.9 倍 (t 2 200 fm/c)。 在 QMD 类 输 运 
模型 计算 中 , 轻 粒 子 产 额 (尤其 是 ?He 和 4He) 通常 低 于 
实验 值 ， 在 我 们 之 前 的 研究 中 时，UrQMD 模型 计算 给 
出 的 3He 的 产 额 大 约 比 实验 值 低 3 倍 。 在 图 2 中 ， 可 以 
看 出 随 着 GEMINI++ 的 加 入 (考虑 衰变 ) 在 一 定 程度 上 
改善 了 对 轻 粒 子 GHe 和 4He) 产 额 的 描述 。 但 是 这 些 轻 
碎片 的 产 额 还 是 低 于 相关 的 实验 值 ; 此 外 ，GEMINI++ 
的 加 入 也 导致 了 自由 核子 产 率 的 高 估 。 为 了 解决 这 些 问 
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AutAu Ey = 59 MeV/nucleon b= 0-2 fm 子 的 集体 流 v, 和 v, MEA TEKE, v, EP OR 
: with GEMINI 度 的 斜率 变 小 及 v, 的 值 增加 。 由 图 2 可 知 ， 考 虑 预 碎 片 
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2 在线 彩 图 ) 束 流 能 量 
缘 碰撞 (b = 0-2 fm) 中 不 同时 刻下 是 
t, *He 和 4He 粒子 的 快 度 分 布 


题 ， 可 以 进一步 在 输 运 模型 中 考虑 粒子 的 自 旋 自 由 度 和 
BFP EAE, EXAT 中 ， 轻 碎片 的 动力 学 产生 
机 制 被 引入 到 玻 尔 兹 曼 类 输 运 模型 中 ; 在 文献 [所 
轻 碎 片 的 直接 产生 道 被 引入 到 QMD 类 输 运 模型 中 。 

些 研 究 表明 轻 碎片 的 动力 学 产生 机 制 对 合理 描述 轻 碎 B 
的 实验 数据 十 分 重要 。 


否 考 虑 衰变 后 p, d, 


3.2 集体 流 


集体 流 作 为 常用 的 核 物质 状态 方程 的 敏感 观测 量 ， 
其 信息 将 对 精确 提取 核 物质 状态 方程 的 信息 有 着 关键 影 
响 M71。 接 下 来 我 们 将 探究 次 级 衰变 对 轻 粒 子 的 集体 
流 分 布 的 影响 及 其 作用 机 制 。 图 3 给 出 了 入 射 能 量 为 
40 MeV/nucleon 时 Au+Au 碰撞 中 《He 粒子 的 直接 流 v, 
(b = 2 一 5.5 fm) 和 椭圆 流 ww (b = 5.5 一 7.5 fm) 的 约 化 快 
度 分 布 ， 演化 停止 时 间 为 += 200 fm/e。 当 考虑 了 预 碎片 
的 退 激 之 后 ， 新 产生 的 “He 粒子 主要 来 源 于 大 质量 碎片 


的 衰变 之 后 ,，“He 粒子 的 产 额 有 着 
布 出 现 了 双 峰 结构 ， 相 
子 集体 流 的 分 布 产生 一 定 的 影响 。 
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得 到 的 结果 更 接近 实验 测量 值 。 因 为 重 离子 们 撞 过 
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(增加 )。 这 是 因为 在 t= 150 fme 附近 ， 核 子 处 于 反应 末 


态 的 低 密 环境 中 ， 此 时 核子 之 
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的 大 质量 碎片 存在 ， 对 了 
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组 的 结果 T2080 (E 
碎片 衰变 对 vi 的 影响 较 弱 
变 对 轻 粒 子 对 集体 流 的 影响 


) 对 比 


AHAT, uo > 0.8) 时 不 同 碰撞 参数 下 得 到 的 轻 粒 子 
化 的 结果 及 与 INDRA 和 FOPI 实验 


。 通 过 对 比 这 两 个 能 量 下 衰 


可 以 发 现在 碰撞 能 量 相对 较 


低 时 预 雁 片 的 退 激 效应 更 加 明显 ， 这 是 因为 在 较 低 的 碰 


撞 能 量 下 较 大 质量 的 碎片 更 容易 形成 ， 而 这 些 碎 片 将 导 
致 轻 结 团 产 额 的 升 高 。 


据 轻 粒子 产 额 、 


围 内 ,只 有 在 低能 区 ,由 于 受到 相互 作用 影响 的 时 间 较 长 、 
产生 的 大 碎 块 质量 较 大 等 因素 的 影响 , 相 较 于 在 200 fm/c 
时 停止 输 运 模拟 接 入 GEMINI++ 的 结果 ， 在 300fm/c 时 
考虑 预 碎片 衰变 的 轻 粒 子 的 集体 流 改变 量 大 约 为 2096; 
在 高 能 区 ， 集 体 流 的 变化 仅 为 0.2%。 

图 5 给 出 了 束 流 能 量 Ej, = 40 ~ 400 MeV/nucleon 
范围 内 集体 流 参数 v 和 vo 的 激发 函数 并 与 实验 值 进 
行 了 对 比 。 在 不 考虑 预 碎片 的 衰变 时 ， 较 低能 量 下 的 vj 
明显 高 于 实验 值 ; 然而 , 一 旦 考虑 衰变 v (v) 的 值 明 
显 地 被 压低 (提高 )， 这 与 图 4 的 结果 相 一 致 。 此 外 ， 所 
研究 能 区 范围 内 ，GEMINI++ 的 加 入 使 得 模拟 计算 结果 
v, 和 vo 更 加 接近 实验 值 。 同 时 ， 随 着 束 流 能 量 的 增加 ， 
预 碎片 的 衰变 对 集体 流 的 影响 逐渐 减 小 ， 因 为 随 着 能 量 
的 升 高 ， 反 应 变 得 更 加 剧烈 ， 大 质量 预 碎片 的 多 重 数 降 
低 ，GEMINI++ 处 理 的 激发 核 也 相应 变 少 。 


33 核 阻 止 本 领 
核 阻止 本 领 通常 可 以 用 来 表征 碰撞 核 的 透明 度 ， 反 


II 


况 , E t = 200 fm/c 之 后 
量 的 改变 相对 较 小 ， 故 而 我 们 将 模拟 


fm/c。 这 里 需要 指出 的 
阶段 所 经 历 的 时 间 是 不 一 样 的 ， 但 在 本 文 研究 的 


体 流 和 核 阻止 随时 间 的 演化 情 
米 续 演 化 所 带 来 的 轻 粒子 观测 
的 时 刻 停 在 1 = 200 
有 然 不 同 能 量 下 反应 的 动力 学 
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(a) 


映 了 碰撞 系统 原始 纵向 束 流 能 量 在 反应 中 转化 为 内 部 自 
由 度 的 程度 。 研 究 表 明 核 阻止 本 领 敏感 于 核 物质 状态 方 
程 和 两 体 散 射 截面 ， 因 此 常 被 用 于 核 物质 状态 方程 的 约 
东 以 及 两 体 散射 截面 的 研究 中 B35 和 。 因 此 ， 研 究 衰 变 
对 核 阻 止 本 领 的 影响 将 会 对 状态 方程 以 及 两 体 散 射 截面 
等 相关 研究 有 着 积极 的 作用 。 图 6 展示 了 在 束 流 能 量 为 
40 MeV/nucleon 的 Au+Au 中 心 碰撞 (b = 0-2 fm) 中 考 
E ( 短 划 线 / 实 线 ) 和 不 考虑 (点 虚线 /点 线 ) 衰变 前 后 ?He 
和 “He 粒子 在 x ( 蓝 色 线 ) füz (绿色 线 ) 方向 的 快 度 分 
布 ， 以 及 相应 的 RE 的 值 。 可 以 看 到 衰变 对 ?He 和 ^He 
粒子 的 快 度 分 布 有 着 明显 的 差异 ， 但 是 相应 的 RR, 值 的 
差别 却 很 小 。 结 合 图 2 的 结果 ， 可 以 看 出 预 碎片 的 衰变 
确实 会 影响 轻 粒 子 的 产 额 分 布 , 但 是 由 于 GEMINI++ 中 
的 衰变 在 预 碎片 质心 系 内 是 各 向 同性 的 ， 因 此 对 于 各 向 
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&|5 (在 线 彩 图 )Au+Au 碰撞 过 程 中 提取 的 集 
和 vao (b) 随 束 流 能 


)。 在 


400 


SATB v, (a) 
量 的 变化 及 INDRA™ A] Fop?! 
没有 显示 误差 条 的 地 方 ， 它 


异性 参量 的 核 阻止 本 领 影响 较 小 。 
与 图 5 相同 ， 图 7 给 出 了 Autau 中 心 碰撞 中 考虑 
衰变 前 后 自由 质子 (a) 和 Z = 1 粒子 (b) 的 核 阻 止 本 
领 RE 和 varti 的 束 流 能 量 激发 函数 。 很 明显 ， 考 虑 衰变 
Ja, Re 和 zart 的 数值 均 被 拉 低 ， 并 且 能 更 好 地 描述 实 
验 数据 。 此 外 ， 误 变 对 核 阻止 本 领 的 影响 随 着 束 流 能 量 
的 增加 而 逐渐 减弱 。 因 为 在 相同 的 碰撞 参数 下 ， 随 着 入 
射 能 量 的 升 高 ， 反 应 程度 加 剧 ， 产 生 的 大 质量 激发 态 预 
肆 片 多 重 数 逐渐 降低 ,故而 预 俯 片 的 退 激 效应 逐渐 减弱 。 
与 产生 的 预 轻 碎 片 和 自由 粒子 不 同 ， 由 大 质量 的 激发 态 
预 碎片 衰变 产生 的 轻 碎 片 和 自由 粒子 的 出 射 是 各 向 同性 
的 ， 这 些 较 大 质量 的 预 碎 片 通常 具有 较 弱 的 核 阻止 本 领 
(RE 和 vartl 较 小 )E5， 因 此 衰变 产生 的 轻 碎片 和 自由 
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CNPC2023: Effects of Sequential Decay on Collective Flows and Nuclear 
Stopping Power in Heavy-ion Collisions at Intermediate Energies 


XIAO Kui! ?, LI Pengcheng!, WANG Yongjia!, LI Qingfeng! + 


(4. School of Science, Huzhou University, Huzhou 313000, Zhejiang, China; 
2. School of Physics and Microelectronics, Zhengzhou University, Zhengzhou 450001, Henan, China; 
3. Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, Gansu, China ) 


Abstract: In heavy ion collisions (HICs), the production of light particles plays an important role in extracting informa- 
tion about the equation of state (EoS) of nuclear matter. Based on the Ultra-relativistic Quantum Molecular Dynamics 
(UrQMD) model, the effect of the sequential decay on the collective flows and the nuclear stopping power of light par- 
ticles in Au+Au collisions at intermediate energies were investigated, with the statistical decay model GEMINI++ is 
used to process the secondary decay of the primary fragments. It is found that due to the memory effect and the daugh- 
ter nuclei produced by decay inherent part of the dynamic information of the parent nucleus, the experimental data can 
be better described by considering the sequential decay. And the influence of the sequential decay on the observables 
weakens with the increase of the collision energy. The results highlight the sequential decay and the production of light 
particles in HICs have an obvious effect on the observables which sensitive to the EoS, and these effects should be 
considered when adopting these observables to extract the information of the EoS. 


Key words: heavy-ion collisions; statistical decay; collective flow; nuclear stopping power 
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